An Exactly Solvable Anisotropic Directed Percolation Model in Three Dimensions 
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We solve exactly a special case of the anisotropic directed bond percolation problem in three- 
dimensions, in which the occupation probability is 1 along two spatial directions, by mapping it to 
a 5-vertex model. We determine the asymptotic shape of the infinite cluster and hence the direction 
dependent critical probability. The exponents characterizing the fluctuations of the boundary of the 
wetted cluster in d-dimensions are related to those of the (d — 2)-dimensional KPZ equation. 



It is well-known that the critical behaviour of a very 
large class of reaction-diffusion systems is describable 
in terms of the critical exponents of directed percola- 
tion (DP) 0. Examples include heterogeneous cataly- 
sis contact processes and epidemic models, self orga- 
nized criticality |j and the driven depinning transition 
j| . Though rather precise estimates of critical parame- 
ters are known from various numerical techniques such 
as series expansions and Monte Carlo simulations ||, an 
exact solution has not been possible so far, even in the 
simple case of 1+1 dimensions. The only analytically 
tractable case known is the anisotropic bond percolation 
on a square lattice, solved by Domany and Kinzel (DK) . 
They addressed the problem in which horizontal and ver- 
tical bonds are present with different concentrations, and 
solved exactly the special case when one of these is set 
equal to 1 Q. In this case, the wetted cluster has no 
holes, and this fact helps in reducing the problem to that 
of a simple random walk in one dimension . 

In this paper, we solve the three dimensional gener- 
alization of DK model. We consider the DP on a three 
dimensional simple cubic lattice with bond concentration 
in the three directions being p x , p y and p z , and solve the 
case p x — p y — 1 , p z arbitrary. We study the properties 
of the infinite cluster and obtain the macroscopic shape 
of the cluster and the scaling form for the fluctuations of 
its boundary . By universality, the behaviour of fluctua- 
tions of the outer boundary of the infinite cluster at large 
length scales in this special soluble case is expected to be 
the same as in the more general case above percolation 
threshold with p x ,p y , and p z arbitrary. 

In addition to the usual critical exponents of the DP 
problem, which are defined in terms of the power-law 
behaviour of various quantities near the critical point, 
there are universal power-law prefactors to the exponen- 
tial decay of correlation functions away from the critical 
point. For example, for p < p c , the probability of a finite 
cluster having s sites varies as s~ e exp(— As) as s — ► oo. 
For p > p c , the probability varies as s~ e exp(—Bs~3~) 
for large s. The exponents 9 and 6' are examples of off- 
critical exponents [^j . We find that the fluctuations of the 
boundary of the infinite cluster in d-dimensions are in the 
same universality class as those of a (d-2)-dimensional in- 



terface moving in a (d-l)-dimensional space. The latter 
are described by the well-known KPZ equation ||]. Thus 
we identify the critical exponents of KPZ equation as 
belonging to the class of off-critical exponents of DP. 

Our solution falls in the general class of disorder solu- 
tions of statistical mechanical models |u| . These often 
show the phenomena of dimensional reduction. Thus a 
3-dimensional problem at disorder point can be thought 
of as a 2-dimensional system evolving in time, and cor- 
relation functions in its steady state would show large 
distance tails characteristic of 2-dimensional systems. In 
our case, the choice p x = p y = 1 allows a further reduc- 
tion of dimension by 1, and we show that the system is 
equivalent to describing the Markovian evolution of a sys- 
tem of hard-core particles on a 1-dimensional ring, where 
the particles are able to move in only one direction. [|ll| . 
Our solution is of interest also as an exact solution of a 
nontrivial three-dimensional statistical mechanics prob- 
lem with positive weights, of which not many are known 

B 

Consider the general DP problem on a simple cubic 
lattice with nearest neighbor bonds, all directed in the 
direction of increasing coordinates. The concentration of 
bonds is p x ,p y and p z along the x-, y- and z- axes. We 
imagine a source of fluid at the origin, which can spread 
along the occupied directed bonds. The sites wetted by 
this fluid form a cluster which we shall call the wetted 
cluster for brevity. One would like to calculate the prob- 
ability, P{R), that a site R belongs to the wetted cluster. 
When p x ,Py and p z are small, the wetted cluster is finite. 
As the probabilities are increased, an infinite connected 
path appears. We shall call the special direction along 
which an infinite path first appears as the preferred direc- 
tion. As the bond concentrations are further increased, 
percolation occurs in a narrow cone centered around the 
preferred direction. Following DK, we look at the angu- 
lar dependence of P{R). Its exponential decay for large 
\R\ defines an angle dependent correlation length, £n, in 
the direction Q. 

The correlation length £n diverges as [p c (f2) — p]~ v , 
where p c (f2) is the f2 dependent critical probability. The 
value of the exponent v depends on the direction f2. We 
show below that in most directions v takes the same 
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value 3/2. is different if f2 is the preferred direction 
(when v = z/jj d , the DP exponent), or when it lies in 
the one of coordinate planes (xy,yz and zx ). For per- 
colation in these planes, clearly all perpendicular bonds 
can be ignored, and the problem reduces to the known 2- 
dimensional case. Consider, for simplicity, the isotropic 
case (p x — p y = p z )- Then the easy direction is (111), 
and along this direction vf 1 1.29. Along (110), (Oil) 
and (101) directions v^ d « 1.733. Other directions in the 
xy,yz and zx planes have v — 2. The x-, y- and z- axes 
have v = 1. 

Consider a coordinate system in which x- and y- axes 
are in the plane of the paper with the z-axis pointing 
out of it (see fig. |l|). Since p x = p y = 1, all the bonds 
aligned along the x-axis and y-axis are present. Hence, 
if a point (x,y, z) is wetted, then so are all the points 
(x',y',z) with x' > x,y' > y. It is easy to see that if 
(x,y,z) is wet, so must all the points directly below it 
in the ^-direction. Therefore, we can define an integer 
height function h(x,y) such that all points (x,y, z) with 
z > h(x, y) are dry while the points with z < h(x, y) are 
wet. The wetted cluster has no holes and can be specified 
completely by its bounding surface, h(x,y). 

We first determine h(x,y), the mean value of h(x,y). 
We do so by mapping the problem to a 5-vertex model 
p3| . Consider the orthogonal projection of the wetted 
cluster on to the xy plane (see fig. 0(a)). It consists of 
paths running from y-axis to the x-axis via steps in the 
right and down directions obeying the constraint that 
paths do not cross each other. The k th path separates the 
sites with h(x,y) < k from those sites with h(x,y) > k. 
We now slide each point in the k th path by (fc, k) (see 
fig. ||(b)), thereby mapping the point {x,y) to the point 



This ensures that two paths will not have common edges. 
Now, every site has either zero or two bonds connecting 
it to it's neighbors. Thus a given configuration of shifted 
paths will be made up of 5 kinds of vertices. It is easy 
to check that the correct weights of each of these ver- 
tices are as shown in fig. pi where q — 1 — p. Under this 
mapping, it is an elementary exercise to show that the 
gradients are related through 

h( = h x /(l + h x + h y ) , 

h v =hy/(l + h x + h y ) , (2) 

where stands for keeping the second coordinate rj 
fixed, and similarly for other partial derivatives. 

Consider a point (£,rj) with £,r] 3> 1. In the 5-vertex 
problem, h( and h n are the mean densities of lines in the 
horizontal and vertical directions respectively. They are, 
however, not independent and knowledge of one deter- 
mines the other. The relation between and h v is a 
local relation and is most easily determined using a dif- 
ferent set of boundary conditions when h( and h v are uni- 
form everywhere. Thus we consider the 5-vertex model 
on an open cylinder of length N in the y direction, and 
infinite in the x direction. 

The 5-vertex model can be solved via the transfer ma- 
trix technique. We transfer from column to column along 
the x-axis. Consider the sector in which there are n lines 
in a column. The ice-rule plus periodic boundary condi- 
tions ensure that there are n lines in each column [ p^[ . 
The transfer matrix, T({y'}; {y}), which is the weight of 
going from configuration {y} to {y'} (see fig. ^ where the 
labelling is self explanatory) is the product of the weights 
of each vertex and can be written as 




FIG. 1. A typical wetted cluster shown up to z = 6 (for 
p=0.3). The cluster extends infinitely in the x and y direc- 
tions. 



it 

T({y'};{y}) = I[(p 1 ~ 
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(C, rj) given by 



( = x + h(x,y), 
T] = y + h(x,y). 



(1) 



(a) (b) 

FIG. 2. (a) The projection of the wetted cluster on to the 
x-y plane. Overlapping edges have been shown separated for 
clarity, (b) The shifted paths. 
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q 1 1 p 1 

FIG. 3. The vertex weights. 

It is easy to verify that Yl{ y n T{{y'}; {y}) = 1, where 
the summation is over all allowed configurations. There- 
fore T({y'}; {y}) is a properly normalized transition 
probability of ({y} -> {y'}). 

One can visualize this process as a system of n hard 
core particles, each hopping only to the right on a ring 
of size N. If the particle on the right is at a distance 
m then the particle can hop up to m — 1 steps during 
one time step, with the probability for k steps being 
P(k\m) = p^-SkOqin-k-i _ -p Qr i ar g e times, the above, 

being a Markov process, will evolve into its steady state. 
The steady state will be one in which all states are equally 
likely. To see this, we note that Yl{ y \ T({y'}; {y}) = 1. 
Therefore T 4 , the transpose of T, is also a stochastic ma- 
trix. Hence it has the state (1, . . . , 1) as the left eigen- 
vector with eigenvalue 1. Taking transpose, (1, . . . , 1)', 
is a right eigenvector with eigenvalue 1. 

Knowing the steady state, we can determine the aver- 
age number of steps, d, of a particle in the y direction for 
each transfer in the x direction. Simple algebra gives 



P(l ~ P) 
p(p + IP) ' 



(4) 



where p = n/N. Making the identification d p = h^jh^ 
and p = h v in eq.(^), we get 



p (l - h{ - h v ) = qh^hn 



(5) 



Eq.(||) can be rewritten in x, y coordinates with the help 
of eq.(pl) to give 



p (l + h x + hy) = qh x h y . 
For large x, y, eq.(ph has the scaling solution 



H x ,v) = -h(Ax,Ay). 



(6) 



(7) 



Putting the ansatz h(x, y) = A(x + y) + B^/xy, we 
find that it satisfies the equation with A = p/q and 
B = 2^p/q. 
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FIG. 4. Two typical arrangements of lines in adjacent col- 
umn (for n = 6). The position of lines, {?/;}, in the upper 
column must interlace with the position, {y'i}, in the lower 
column. 



Thus the asymptotic shape of the surface is given by 
h{x,y,p) = (px+py + 2yfpxy)/q, for x, y>l. (8) 

The expression is symmetric in x and y. If y = 0, 
then h = px/q as expected. For given p, and given 
the azimuthal angle, <p, there is a critical polar an- 
gle 9 C such that all points with 8 > 8 C are wetted. 
The critical angle is determined through the relation 
cot(0 c ) = h/^/x 2 + y 2 . Eq.(||) can then be inverted to 
derive the direction dependent critical probability, the 
smallest probability for which (9, <j>) is wetted, as 



V(q + !)(/?+!)- VoTp" 
a + (3 + 1 



(9) 



where a = tan(0) cos(</>) and j3 = tan(0) sin(</>). 

The problem can also be viewed as a 1-dimensional 
cellular automaton evolving in time. Consider the coor- 
dinate transformation t = x + y,u = x~y. Then from 
the definition of the model it is easy to check that for 
t > 0, h(u,t) follows the rule, 

h(u, t + 1) = max[h(u - 1, t), h(u + 1, t)] + Tj(u, t), (10) 

where r](u,t)'s are i.i.d. random variables taking non 
negative integer values with Prob(r) = k) = (1 — p)p k . 
The rules of evolution are local and we expect the fluc- 
tuations to be governed by the KPZ equation Q 

We now look at the correlations in the system. There 
are two correlation lengths of interest. £n(f2) has been 
defined through the exponentially decaying probability 
of wetting. £ii(fi) diverges as [p c (f2) — p]~ v ■ We denote 
by £j_ the correlation length along equal height contours 
but along the wetted surface. We can study the fluctua- 
tions by looking at the scaling properties of the 5-vertex 
model. The transfer matrix can be diagonalised via the 
Bethe ansatz From this analysis it is known that 

the dynamic exponent z — 3/2, and the roughening ex- 
ponent x = 1/2 fl(| • These exponents are the well known 
exponents of the KPZ equation. 

The height fluctuations at a point are linearly related 
to the fluctuations in the trajectories in the 5-vertex 
problem. In the KPZ problem, it is known that the mag- 
nitude of the fluctuations of height varies as /i 1 / 3 . Thus 
we get y/ ((Sh) 2 ) ~ h 1 ^ 3 . £|| is the correlation length 

at an angle 6 C + 59. Clearly ^69 - h 1 / 3 - or 

f|| - (S9)~ 3/2 . From eq.(§), we note that, for p ^ 0, 
59 ~ 5p. Thus we get v — 3/2, to be compared with the 
value v = 2 for the two-dimensional DK. The behaviour 
of £x can be obtained from the relation £_l 



we get that £j_ 

height contour from the origin. 



2/3 



£j /z . Thus 

, where arii is the distance of the 
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Other 3-dimensional lattices can be treated similarly. 
We consider lattices made up of two-dimensional layers 
stacked on top of each other with vertical bonds con- 
necting a site to the one vertically above it. If each 
layer is a triangular lattice, then it may be described as a 
square lattice with one set of diagonal bonds. However, 
the diagonal bonds provide no additional connections if 
Px = Py = 1) and can be neglected. Thus the problem re- 
duces to the square lattice. If each layer is a honeycomb 
lattice, there are two types of sites: those having only 2 
or 3 outgoing bonds. If we integrate over the former, the 
problem again reduces to a simple cubic lattice, with a 
renormalized probability p z . 

In higher dimensions, d > 3, the exponents are re- 
lated to the higher dimensional KPZ equation. The 
d-dimensional DK model is equivalent to a (d — 2)- 
dimensional interface growing in time. As in the case of 
depinning transition of a driven tilted interface |i~7] l , the 
height-height correlation function in e?-dimcnsions will 
have the scaling form 

(5h(x)Sh(x')) = I*,, -x\\^F f 1 ^— -M- \ , 

(11) 

where Xn is measured along the radial direction, x± is 
measured along the equal height contour and Zd and Xd 
are KPZ exponents of a growing <i-dimensional interface. 

In summary, we have been able to exactly solve the 
three-dimensional anisotropic directed percolation prob- 
lem for the special case p x — p y = 1. We used the 
fact that in this case the wetted cluster has no holes 
within. This reduces the problem to studying the two- 
dimensional surface of the cluster. The weights of differ- 
ent surfaces were shown to be same as those of different 
configurations of a 5 vertex model. From the already 
known exact solution of the latter, we obtained an exact 
expression for the average height profile h(x, y) and de- 
termined the asymptotic correlations of the surface fluc- 
tuations of the wetted cluster for large separations. 

We thank Mustansir Barma and Satya Majumdar for 
critical reading of the manuscript. 



[1] P. Grassberger, J. Stat. Phys. 79, 13 (1995). 



[2] R. M. Ziff, E. Gulari, Y. Barshad, Phys. Rev. Lett. 56, 
2553 (1986). 

[3] S. Maslov, Y. C. Zhang, Physica A. 223, 1 (1996); B. 

Tadic and D. Dhar, Phys. Rev. Lett. 79, 1519 (1997). 
[4] L. H. Tang and H. Leschhorn, Phys. Rev. A. 45, R8309 

(1992); S. V. Buldyrev, A. L. Barabasi, F. Caserta, S. 

Havlin, H. E. Stanley and T. Vicsek, Phys. Rev. A. 45, 

R8313 (1992). 

[5] I. Jensen, J. Phys. A. 29, 7013 (1996); J. W. Essam, A. J. 
Guttmann and K. De'Bell, J. Phys. A. 21, 3815 (1988); 
R. J. Baxter and A. J. Guttmann, J. Phys. A. 21, 3193 
(1988); P. Grassberger and Y. C. Zhang, Physica A. 224, 
169 (1996). 

[6] E. Domany and W. Kinzel, Phys. Rev. Lett. 47, 5 (1981). 
[7] F. Y. Wu and H. E. Stanley, Phys. Rev. Lett. 48, 775 
(1982). 

[8] A recent discussion of the exponent 9 can be found in S. 
Lai and M. E. Fisher, J. Chem. Phys. 103, 8144 (1995). 
For undirected percolation, within droplet model, the 
value of 9' is known in all dimensions (see T. C. Lubensky 
and A. J. McKane, J. Phys. A. 14, L157 (1981)). 
[9] M. Kardar, G. Parisi and Y. C. Zhang, Phys. Rev. Lett. 
56, 889 (1986); T. Halpin-Healy and Y. C. Zhang, Phys. 
Reports, 254, 216 (1995). 

[10] P. Rujan, J. Stat. Phys. 49, 139 (1987); A. Georges and 
P. Le Doussal, J. Stat. Phys. 54, 1011 (1989). 

[11] In the case of off-critical directed animal problem, the 
dimensional reduction by 1 is discussed in J. L. Cardy, J. 
Phys. A. 15, L593 (1982). 

[12] For other soluble models in 3-dimensions, see V. V. 
Bazhanov and R. J. Baxter, J. Stat. Phys. 69, 453 (1992); 
H. Y. Huang, V. Popkov and F. Y. Wu, Phys. Rev. Lett. 
78, 409 (1997); V. B. Priezzhev, J. Stat. Phys. 26, 817 
(1981); J. F. Nagle, C. S. O. Yokoi and S. M. Bhattachar- 
jee in Phase Transitions and Critical Phenomena Vol 13, 
cd. C. Domb and J. L. Lebowitz (Academic Press, Lon- 
don ,1990); M. Prahofer and H. Spohn, J. Stat. Phys. 88, 
999 (1997). 

[13] This is a special case of the 6-vertex model for which 
the exact solution is known. For a review of the 6-vertex 
model, see E. H. Lieb and F. Y. Wu in Phase Transitions 
and Critical Phenomena Vol.1, ed. C. Domb and M. S. 
Green (Academic Press, London and NY, 1972), pp 331- 
490. References to more recent work may be found in H. 
Y. Huang, F. Y. Wu, H. Kunz and D. Kim, Physica A. 
228, 1 (1996). 

[14] p 363 in Lieb and Wu op. cit. 

[15] pp 374-375 in Lieb and Wu op. cit. 

[16] D. Dhar, Phase Transitions. 9, 51 (1987); L. L. Gwa and 

H. Spohn, Phys. Rev. A. 46, 844 (1992). 
[17] L. H. Tang, M. Kardar and D. Dhar, Phys. Rev. Lett. 

74, 920 (1995). 



4 



